Subgap conductivity in SIN-junctions of high barrier transparency 
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We investigate the current-voltage characteristics of high-transparency superconductor-insulator- 
normal metal (SIN) junctions with the specific tunnel resistance p < 30 Q x pm 2 . The junctions 
were fabricated from different superconducting and normal conducting materials, including Nb, Al, 
AuPd and Cu. The subgap leakage currents were found to be appreciably larger than those given by 
the standard tunnelling model. We explain our results using the model of two-electron tunnelling 
in the coherent diffusive transport regime. We demonstrate that even in the high-transparency 
SIN-junctions, a noticeable reduction of the subgap current can be achieved by splitting a junction 
into several submicron sub-junctions. These structures can be used as nonlinear low-noise shunts in 
Rapid-Single-Flux-Quantum (RSFQ) circuitry for controlling Josephson qubits. 



INTRODUCTION 

In the past decades, the mechanisms of the conductiv- 
ity through the superconductor-insulator-normal metal 
(SIN) tunnel interface have been extensively studied. 
Due to the superconducting energy gap A in the quasi- 
particle density of states, the current through such an 
interface is generally suppressed at low temperatures, 
fceT -C A, and subgap bias voltages, V < V g = A/e, giv- 
ing rise to a low- voltage nonlinearity in the TV-curve [l| . 
Recently this, the most common property of SIN-contacts 
has been found to be rather useful, e.g., for reducing gen- 
eration of quasiparticles by building the current shunts 
into the Josephson-junction qubits [2j or for avoiding an 
extra decoherence Q in all- Josephson RSFQ-qubit inte- 
grated systems 0. 

In particular, the thermal noise produced by the resis- 
tively shunted Josephson junctions in a standard RSFQ- 
circuitry, brings into being an additional source of qubit 
decoherence 3- It was shown however that this effect 
can be minimized, by using the nonlinear damping that 
is provided by an SIN-junction of low asymptotic resis- 
tance, i?N, and a high value of the nonlinearity parameter 
i] = [G(0).Rn] _1 , where G(0) is the zero-bias conductiv- 
ity of the shunt. The dynamics of SIN-shunted Josephson 
junctions has recently been analyzed in detail on the ba- 
sis of a standard tunnelling model of the SIN-junction 
(see Ref. @). The feasibility was demonstrated to en- 
sure both a high damping at the characteristic Joseph- 
son frequency u>c = 2eVc/h, where Vq — IcRn is the 
characteristic voltage and Iq the critical current, and, 
due to the high impedance of the SIN-junction at subgap 
voltages, a sufficiently low noise at qubit frequencies of 
typically up to - 10 - 30 GHz < A/h. 

The first experiments on SIN-shunted Nb SIS Joseph- 
son junctions [3j have proven the feasibility of the over- 
damped regime at T = 4.2 K. At the same time, how- 
ever, an unexpectedly weak subgap nonlinearity of the 
SIN-junctions was observed. In particular, the value of 
?/ « 6.2, measured at T = 1.4 K in an SIN-junction 
with a specific tunnel resistance of p ss 135 O x /xm 2 , 



was much smaller than the values r\ > 100 which are 
typical of opaque barriers (see, e.g., Refs. 0, and 
also far below the simple estimates made for the tun- 
nelling model, yielding a large subgap suppression fac- 
tor cx exp(— A//cbT) [l|. In view of the previous data, 
it would be an even more challenging task to achieve a 
high nonlinearity in SIN-junctions of very low specific 
resistance, e.g., p < 30 57 x pirn 2 for the SIN-junctions 
based on superconducting Al [3j, as required for RSFQ 
operation at qubit temperatures T < 50 mK. 

In this paper, we address subgap conductance phenom- 
ena in SIN-junctions of very high transparency. In par- 
ticular we found that the degree of nonlinearity of the 
measured 7 V-curves is consistent with predictions of the 
two-electron tunnelling model which has been developed 
in a ballistic regime by Blonder, Tinkham and Klapwijk 
(BTK) 7'] and modified by Hekking and Nazarov |8|,|9j to 
account also for interference contributions due to coher- 
ent electron diffusion. In earlier experiments, the spatial 
coherence was clearly demonstrated by Pothier et al. 
in the form of a periodical dependence of the low-bias 
conductance of the so-called NS-QUID interferometer on 
the magnetic flux. A spatially coherent enhancement of 
the two-electron tunnelling amplitude makes the junc- 
tion conductance topology-dependent. This effect was 
observed as zero-bias conductance anomaly in a thin- 
film planar NS-interface, which was compared in Ref. fll| 
with the BTK ballistic transport through an edge-type 
NS-contact. In our experiment, we compare the val- 
ues of the nonlinearity parameter r\ for different junction 
topologies, varying the junction sizes and/or the width 
of the adjacent diffusive electrode. We show that a suffi- 
ciently strong nonlinearity can be obtained in submicron 
SIN-contacts with even high-transparency barriers, i.e. 
in the case of interest for RSFQ-qubit applications. 



DESIGN OF SAMPLES 

Two techniques were applied for fabricating the sam- 
ples. The large reference junctions "NbAuPd" and 
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FIG. 1: SEM-micrographs of the sample 1L1F (a) and blow- 
up of the junctions 1L2F (b), 1L4F (c), 2L2F (d), and 
lLlFmin (e), of different geometries. The total tunneling 
area is the same for each sample, producing similar asymp- 
totic resistances (except for the "2Lx" -samples [(d) in this 
Figure], as measured in a 4-point configuration. 



"NbAl" (see Table QJ were made in a standard sandwich 
technique which is usually applied for RSFQ-devices (see, 
e.g., Refs. |l2|,|l3j])i on the basis of a Nb trilayer, pat- 
terned with the help of optical lithography. Both samples 
had a sophisticated multilayered structure which made a 
particular analysis of the SIN-junction component diffi- 
cult. For example, the junction "NbAl" consisted of the 
metal stack with a pronounced proximity effect between 
the AuPd and the top Nb layers, and, possibly, between 
the bottom Nb and the covering Al, separated by a thin 
oxide barrier only. These Nb-based samples were charac- 
terized in a He -cryostat at temperatures down to 1.4 K. 
The data for these two samples, shown in Table HJ repre- 
sent our rough estimations. 

For the sake of simplicity and flexibility of the design, 
all other junctions were fabricated by means of Dolan's 
shadow evaporation process |l4j . The complete struc- 
tures were deposited in-situ from the e-gun through the 



single e-beam-patterned PMMA/Ge/Copolymer mask 
with free-hanging bridges. The junctions were formed 
on the overlapping areas of the first and the subsequent 
shifted layers, deposited through the same openings in 
the mask but at different tilting angles of the wafer in re- 
spect to the incident material flow. The TV-curves were 
measured at T < 0.05 K in a dilution refrigerator. 

The bottom layer of all shadow-evaporated junctions 
was made of aluminum, with a thickness of g?b = 30 nm, 
and oxidized in a very mild regime (see below) to form 
a thin tunnel barrier. To be able to put several dif- 
ferent structures on the same chip, we used besides a 
two-layer also a three-layer design, with Cu-, Cu/Cu- or 
Al/Cu-layers building the top electrodes of the junctions. 
The junctions "SINTOP" and "SINBOT" were simple 
two-layer structures with the tunnel areas differing by 
the factor of two, both oxidized for 5 min at an oxy- 
gen pressure of Po 2 = 1 Pa, which resulted in a specific 
barrier resistance p — 27 f2 x /im 2 . A somewhat more 
complicated layout was used for the junctions " 1L1F" to 
" 2L4F" which were fabricated on the next wafer. The top 
junction electrode was made either of a single layer [the 
top Cu film in the "lLxx"- junctions, see, e.g., Fig.^a- 
c),(e)] or of a double-layer composition [a composition of 
the middle and the top Cu films in the "2Lxx" -junctions, 
see, e.g., Fig. ^d)], evaporated in two consecutive steps, 
with an additional tilting of the wafer in between. Ac- 
cordingly, two different thicknesses of the top electrodes, 
c?t = 30 nm and 60 nm, were obtained. This wafer was 
oxidized for 5 min at lower pressure, Po 2 = 0.1 Pa, in 
order to obtain an even more transparent barrier. Al- 
though all junctions had the same Al bottom layer with 
the same oxidation step, two different values of p were 



(a) SIS'- junction 
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FIG. 2: Cross-section sketch of the samples SISN (a) and 
SISNSIN (b). In the sample SISNSIN, the junction was a 
composition of an SIN and an SIS' sub-junction, both having 
the same Al-I bottom layer (see text). 
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found: p = 1517 x pm 2 for the " lLxx" -junctions, and 
p = 5 f2 x pm 2 for the "2Lxx" -junctions. A similar effect 
was observed for the samples "SISN" and "SISNSIN" . In 
our opinion, this variation is specific to low-dose oxida- 
tion processes and arises due to a prolonged oxidation of 
the Al by adsorbed oxygen in the areas which are not 
being covered by the middle layer of Cu. 

The junctions "1L2F" to "2L4F", "SISN" and 
"SISNSIN" were designed in a "split-finger" layout [see 
Figs.HJb-d)] by subdividing the junction into N — 2 or 4 
narrower sub-junctions with a partial width W each, as 
listed in Table [Q The neighboring fingers of Cu formed, 
together with the superconducting Al sections, the closed 
loops of the area ~ 1 x 1 pm 2 . In contrast to the other 
samples, the normal metal electrodes in the junctions 
" lLlFmin" and " lL2Fmin" were made in a special "elec- 
tron confinement" shape [Fig. Qfe)]: The width of the 
adjacent wire was reduced further to 1/4 of W while the 
tunnel area was kept identical to that in the samples 
"IMF" and "1L2F", respectively. 

The top electrode of junctions in the samples "SISN" 
and "SISNSIN" , both fabricated on the same wafer, con- 
sisted of two different metals, Al and Cu, which formed 
slightly shifted shadows in such a way that, in the sam- 
ple " SISN" , the junction was designed to be effectively of 
type "superconductor-insulator-superconductor" with a 
proximity-reduced gap (SIS'). In the sample "SISNSIN", 
the junction of the same total area was made up of such 
an SIS' to the half with an Al/oxide/Cu SIN-contact (see 
the cross-sections in Fig.|2J. 



RESULTS AND DISCUSSION 

Typical TV-characteristics of high-transparency SIN- 
junctions are shown in Fig. [3] As expected, significant 
subgap leakage currents were observed in most of the 
SIN-junctions, even at the temperature T = 50 mK <C 
A/fcs, resulting in a strong smearing of the supercon- 
ducting gap corner and in rather moderate values of 
i] ~1-10 (see Table [J. The appreciable excess subgap 
current at low temperatures is commonly attributed to 
the resonant two-electron transport channel due to An- 
dreev reflection 15]. To interpret our data, we first make 
estimations within the simplified ballistic BTK model, 
which put aside the electron scattering effects in the elec- 
trodes. Using Eq.(17) of Q, it is straightforward to de- 
rive simple expressions for the subgap conductance at 
T = 0: 



RnG{V) 



2A 2 (l + Z 2 ) 



(eV) 2 + [A 2 - (eU) 2 ] [1 + 2Z 2 } 2 
and the BTK nonlinearity parameter: 

(1 + 2Z 2 ) 2 
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FIG. 3: Comparative JVC-plot for the samples "1L1F", 
" 1L2F" , and " 1L4F" . Whereas the asymptotic slopes of all 
three samples are similar, the zero-bias resistance of the sam- 
ple "1L4F", which is the one with the narrowest wires and 
junctions, exceeds Rj of the other samples noticeably, indi- 
cating the geometric effect for the two-electron conductance. 
The small kicks on both, positive and negative branches of the 
IV-curves around V = ±A/e arise due to self- heating effects 
in the junctions (see, e.g., Ref. (l7jl h 



whereby the dimensionless barrier strength Z |7j char- 
acterizes the transparency of the NS interface, Z = 
kpH/2eF, £f and k F are the Fermi energy and wavenum- 
ber, respectively and H enters the phenomenological re- 
pulsive barrier potential H5(x). For a thin tunnel barrier 
of the simplest, rectangular shape with a height U > e F 
and a thickness a ~ 2nkp 1 = Ap, which is the Fermi 
wavelength, Z can be expressed as 



k F U 
2e F 



sinh (an) 



U -e F 



cf 



(3) 



by matching the well-known plane- wave transmission co- 
efficient (see, e.g., 16]) with that of BTK, C = (1 + 
Z 2 ) -1 , which is provided in Ref. for electrons in 
the normal state. The straightforward evaluation us- 
ing Eqs. J5J and Q and the parameters of the idealized 
Al/A10 x -barrier, U - e F - 2eV (see, e.g., Refs. UlEl) 
and a w 5 A (see the discussion below), as well as 
ep = 11.4eV and fc F = 1.75 x 10 10 m _1 , which are typi- 
cal for Al, results in the values Z « 27 and 7/btk ~ 1500 
and predicts the suppression of the subgap conductivity 
much stronger than found in our experiment. 

Important corrections can be made to the calculations 
if we take a note of the fact that the barriers under consid- 
eration are as thin as a single monolayer of Al oxide. The 
value of ao accepted above presents an average barrier 
thickness of the plane junction with the specific capac- 
itance c « 75 fF /pm 2 (derived from the Coulomb offset 
voltage of the similarly oxidized single-electron transis- 
tors [2(J, cf. also H3), and the dielectric constant e r w 4 
[l9) . For such thin oxide layers, it is plausible to ex- 
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TABLE I: Parameters of the junctions 
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NbAuPd 


Nb /AlOx/AuPd 
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i so i nn 

lOU; lUv 
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380 
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SINTOP 


Al/A10 x /Cu 


0.5x2 


30.50 


0.03 


27 


27 


140 


5.1 


SINBOT 


Al/AlOx/Cu 


0.5x1 


30,50 


0.03 


54 


27 


360 


6.6 


1L1F 


Al/A10 x /Cu 


0.24x2 


30,30 


0.05 


29 
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90 
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lLlFmin" 


Al/A10 x /Cu 


0.24x2 


30, 30 


0.05 
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15 


64 
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1L2F 


Al/A10 x /Cu 


2x(0.24xl) 


30,30 


0.05 


28 
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lL2Fmin a 


Al/A!O x /Cu 


2x(0.24xl) 


30,30 


0.05 


42 


15 
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2.6 


1L4F 


Al/A10 x /Cu 


4x(0.24x0.5) 


30,30 


0.05 


31 


15 
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5.7 


2L2F 


Al/A10 x /Cu 


2x(0.24xl) 


30,60 


0.05 
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Al/A!O x /Cu 


4x(0.24x0.5) 
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3.4 


SISN 


Al/AlOx/Al/Cu 


2x(0.25x0.5) 


30,30,60 
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44 


11 
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SISNSIN 


Al/A10 x /Al/Cu+Al/A10 x /Cu 


2x(0.12/0.12 ! 'x0.5) 


30,30,60 


0.026 


63 


11/22 6 


3.3k 


52 


NbAl 


Nb /A10 x /Al/A10 x /AuPd /Nb 


100 


90,100,100,500 
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1.4 


150 


1.4k 


10 3 



"In these samples, the normal electrode is of "electron- 
confinement" shape (see the text). 
6 For SISN- and SIN-partial junctions, respectively. 



pect considerable local deviations of the relative barrier 
thickness, even for the cleanest tunnel interface. 

Possible models for the fine structure of very thin bar- 
riers have recently been discussed for junctions of type 
Nb/A10 x /Nb with current densities in a wide range, 
J c - 0.1 - 20kA/cm 2 [H |U 0. In Refs. [U l^, 
the barrier is supposed to consist of different areas with 
a certain number of oxide monolayers, including a sta- 
tistically large amount of local pinholes of lateral size 
~ 1 nm, with no single monolayer of oxide at all, domi- 
nating in the current of the whole junction. In contrast 
to this, the model of a uniform, continuously growing 
barrier was suggested by Maezawa et al. [3| as a more 
realistic, predicating upon that oxide is formed on top of 
a granular Al film with a locally different surface orien- 
tation of the crystal planes. Our argument is that small, 
but finite values, Z > 1, corresponding to our experimen- 
tal value r] > 2 (see Fig. 6 of Ref. rather relate to the 
picture of a thin, but finite tunnel barrier, at least on the 
mesoscopic scale ~ L, W . Moreover, on the microscopic 
scale ~ Af, the ballistic transport model accounting for 
the pinhole effects is not adequate; the relevant approach 
here relates to the electron diffraction, which makes the 
tunnelling picture non-local [2^ |. 

To build a quantitative model, we assume a uniform 
distribution of the barrier thickness around the mean 
value ao, i.e. <io(l — 8) < a < oo(l + 8), with a con- 
stant probability density u> = 1/25, and consider an in- 
dependent electron transport through the areas within 
the junction with different local barrier strengths. The 
effective zero-bias nonlinearity can then be found as the 
ratio rj av of the average asymptotic (Rjj 1 ) and the zero- 
bias conductivities (G(0)}: 

„ (ft+gq (4) 

" av "(G(o))-2^ (1 + 2Z2r2 y (4) 

best matching the experimental range of r\ ~ 2-7 in Al- 



based SIN junctions for 5 ~ 0.8-0.9, i.e., for the realistic 
barriers with a thickness of up to two oxide monolayers. 

So far we have been able to obtain a reasonable, though 
rough, estimation of the zero-bias nonlinearity of SIN 
junctions within the scope of the ballistic BTK model. 
A more detailed insight can be gained, by comparing the 
subgap transport properties in SIN junctions of the same 
transparency but of different size. For example, in Fig. [3] 
the voltage-current characteristics of the single-junction 
structure "1L1F" [Fig. [Ja)] is compared to those of the 
split-finger interferometer structures " 1L2F" and " 1L4F" 
[Fig.^b,c)] . All three structures have the same total tun- 
nel area and, therefore, similar asymptotic resistances, 
but they differ with regard to the number N and the 
width W of the partial junctions. In none of the split- 
finger interferometer structures, a periodic dependency of 
the TV-curves on the magnetic field was observed, which 
indicates the absence of inter-junction coherence effects 
such as those reported earlier for different structure di- 
mensions by Pothier et al. 10]. This property of our 
structures is obviously due to the normal-metal sections 
of each loop being, in total, longer than the interference 
length in a diffusive conductor Lt ~ 1 A™ 8] and, when 
used in practice, it makes damping on the base of the 
split-finger SIN-junction flux insensitive. Independent 
tunnelling in the sub-junctions enables a parallelization 
of the SIN-contacts, when a low shunting resistance is 
required. 

For the larger junctions "1L1F" and "1L2F", with 
W ^ 1 /xm, the value of rj = 3.1 does not vary with the 
junction width. On the contrary, for the parallel con- 
nection of the four smaller, 0.5 /zm-wide sub-junctions in 
sample "1L4F", the value of r) = 5.7 is nearly twice as 
large. The same tendency is to be seen for the samples 
"2L2F" and "2L4F" with slightly different parameters 
(see Table HJ . The observed dimensional effect can be 
interpreted following the model of Hekking and Nazarov 
0, which accounts for spatially coherent contributions 



5 



0.6 
0.4 



5~ 0.2 
o 0.0 



.p 

> -0.2 



-0.4 
-0.6 





T = 50mK; B = y 


— 1L1F: n = 3.1 




- - - 1L1Fmin: ri =2 

- 




*~*> r 

/Is j 


\ 
\\ 


1 

>«JV<^* // 




" \V i/ 









- 120 g 



160 



80 £ 



40 



-10 



10 



Current l bias (uA) 



FIG. 4: Current-voltage and differential resistance plots, com- 
parative for the samples "1L1F", with the full- width normal 
wire, and "lLlFmin" with the narrowed normal wire ("con- 
finement" shape). 
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the samples "SISN" and "SISNSIN". A small magnetic field, 
applied perpendicularly to the sample, suppresses the critical 
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to the Andreev conductance due to two-electron diffu- 
sion. According to this model, the coherence factors, and 
thus the subgap conductance, are expected to increase 
with the junction size W until saturation occurred at 
W ~ Lt- For small junctions of size W < Lt, this effect 
can be evaluated using the relation Eq. (12) of Ref. 
which was derived for a normal wire in a tunnel contact 
to a superconducting sheet, rewritten in the form: 

r, coh = [EnG(O)]" 1 cc , (5) 

where ctd is the specific conductivity of the normal elec- 
trode. In our experiment, the validity of Eq. JSJ was 
confirmed by a simple comparison of the nonlinearity pa- 
rameter values in the junctions "1L2F" (ry = 3.1), and 
"2L2F" with 77 = 2 [the layout shown in Fig. [Jd)] with 
a single and a double layer of Cu in the top electrodes, 
correspondingly. The ~ 1.5 times larger value of r\ for 
" 1L2F" arises as the ratio 3 of the specific barrier resis- 
tances, multiplied by the ratio 1/2 of the thicknesses of 
the top electrodes. 

Our data also clearly demonstrate the effect of the 
"confinement" shape for the normal wire. As seen in 
Fig. 01 the zero bias differential resistance of such a sam- 
ple "lLlFmin" [shown in Fig. ^e)] is ~ 1.5 times lower 
than that of the "plain" junction " 1L1F" , even despite its 
higher value of i?N- This tendency also shows in the sam- 
ples " 1L2F" and " lL2Fmin" . From this, we can conclude 
that those configurations of the top electrode where the 
diffusion volume for the back-scattered electrons is re- 
stricted to the closer vicinity at the tunnel barrier result 
in a larger Andreev conductance, which is consistent with 
the predictions of the model Ref. 0, . 

A considerably higher nonlinearity of the TV-curves, 
shown in Fig. O was observed in the junctions "SISN" 
and " SISNSIN" , where the top electrode was made of a 



smaller-gap superconductor, A' <C A, consisting of su- 
perconducting Al, covered by a thicker layer of the nor- 
mal conductor (Cu), see Fig. [2 Due to the second small 
gap in the density of states, the two-electron conductance 
is suppressed at small voltages eV < A' (see, e.g., Q), 
while the multiple Andreev reflections are expected to 
decrease at small bias as the processes of the higher or- 
der in tunnel transparency [22|. In experiment, a strong 
current reduction was found at voltages below ~ A/e 
(cf. also [25| ), resulting in larger values of i] = 43 and 52 
in the samples "SISN" and "SISNSIN", respectively. A 
higher ratio in the latter sample can be explained by the 
lower barrier transparency in the SIN partial junction, 
see Eq. [5] A small critical supercurrent I c < 300 nA, 
otherwise observed, could be easily suppressed below the 
level of fluctuations by a applying perpendicular mag- 
netic field of B ~ 5 mT. 

A much larger value 77 = 10 3 was measured in the sam- 
ple "NbAl" , at the temperature T = 1.45 K which slightly 
exceeds the superconducting transition point for Al. The 
layer of Al was originally considered in this structure to 
form a normal-metal electrode whereas the " S" electrode 
was a proximity superconductor built of an AuPd/Nb 
bilayer. However, much higher subgap nonlinearity than 
in the other SIN-junctions indicates, in our opinion, that 
the Al electrode has a spectrum of states which differs 
from that of a normal metal. For example, a small gap 
could be induced in Al by the bottom film of Nb through 
a very transparent separating tunnel barrier, making the 
junction effectively an SIS' type. We believe that the 
effect of nonlinear damping in SIS structures deserves a 
separate investigation. 



6 



CONCLUSIONS 

We investigated the nonlinearity of the TV- 
characteristics of SIN-junctions of high transparency, 
fabricated with different topologies and composed of 
various materials. Our data are consistent with the pre- 
dictions of the two-electron tunnelling model, accounting 
for the interference contributions in the diffusive scatter- 
ing regime. A clear enhancement of the nonlinearity was 
observed in the submicron junctions with planar dimen- 
sions smaller than the length of interference Lt ~ 1 /^rn. 
We show that in a split-finger junction layout, tunnelling 
in partial junctions is independent and flux-insensitive. 
This property allows a parallelization of these nonlinear 
and, therefore, low-noise submicron SIN-junctions. A 
significant increase in the subgap resistance - compared 
to the SIN-junctions of the similar transparency - was 
observed in proximity-type SIS'-junctions. This could 
be the subject of a special study. 

The values obtained for the zero-bias and the asymp- 
totic junction resistances allow the SIN-junctions to be 
implemented in RSFQ-Qubit circuits as low-noise shunts. 
For example, our junctions with, typically, i?N ~ 10 — 
30 f2 could provide a sufficient damping for submicron 
Josephson junctions with Iq ~ 10 fiA. On the other 
hand, sufficiently long decoherence times, e.g., > 10 [is, 
as estimated in Ref. |2(j for a tunable flux qubit, coupled 
to the RSFQ control circuit, can be achieved due to the 
high zero-bias resistances, [G(O)] 1 ~ 50 — 100 il, of most 
of the measured SIN-junctions. 
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